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It is investigated whether suitable generalizations of the field equations of current field theories to equa¬ 
tions of higher order may be of help in eliminating the divergent features of the present theory. It turns 
out to be difficult, if feasible, to reconcile in this way the requirements of convergence, of positive definite¬ 
ness of the free field energy, and of a strictly causal behavior of the state vector of a physical system. Progress 
may perhaps be made by relinquishing the condition of unlimited localizability of any space time event. 


I. INTRODUCTION 

C ONSIDERABLE progress has been made during 
the last two years in our understanding of the 
scope and limitations of the present picture of systems 
of interacting particles and fields. Thus, noting that a 
theoretical subdivision of the experimental mass and 
experimental charge of the electron so far seems 
unwarranted from any observational point of view, a 
“renormalization program” has been worked out for 
the interaction of electrons and electromagnetic fields. 
This scheme has had great success in making such 
effects as the Lamb-Retherford shift of lines in the 
hydrogen atom and subtle anomalies in the magnetic 
moment of the electron amenable to theoretical inter¬ 
pretation. 1 Nor is the applicability of this evaluation 
technique limited to any power in a development in 
terms of the fine structure constant. 2 

Much less satisfactory, however, has been the out¬ 
come of investigations on the application of renormal¬ 
ization ideas to other systems of particles and fields. 3 
In the domain of nuclear phenomena in particular one 
still encounters many difficulties. Not only do calcula¬ 
tions of certain so-called reactive effects still yield 
divergent results, 4 but even if finite, one often finds 
results which are quantitatively inconsistent with 
experimental data. A striking instance is, for example, 
the ratio of the “extra” magnetic moment of the proton 
to the difference between this moment and the magnetic 
moment of the neutron. 5 It is true that such results 

1 S. Tomonaga et al ., Prog. Theor. Phys. 1, 27 (1946); 2, 101 
(1947); 4, 47, 121 (1949). J. Schwinger, Phys. Rev. 74, 1439 
(1948); 75, 651 (1949); 76, 790 (1949). R. Feynman, Phys. Rev. 
74, 1430 (1948); 76, 749, 769 (1949). 

2 F. J. Dyson, Phys. Rev. 75, 486, 1736 (1949). 

3 See P. T. Matthews, Phys. Rev. 76, 1254 (1949); Phil. Mag. 
41, 185 (1950); also D. Feldman, Phys. Rev. 76, 1369 (1949) for 
electromagnetic properties of vector mesons. 

4 As, e.g., the nucleon magnetic moment calculated on the basis 
of the vector meson theory [see K. Case, Phys. Rev. 75, 1440 
(1949)]. In our present state of limited knowledge it seems 
premature to conclude from the occurrence of such infinities to 
the non-existence of certain particles or certain interactions 
concerned. 

5 K. Case, Phys. Rev. 76, 1 (1949); J. Luttinger, Helv. Phys. 
Acta 21, 483 (1948); M. Slotnick and W. Heitler, Phys. Rev. 75. 
1645 (1949); S. Borowitz and W, Kohn, Phys. Rev. 76, 818 (1949), 


cannot be considered definitive in view of the dubious¬ 
ness of the power series approximation underlying these 
theoretical derivations. Yet the generally unsatisfactory 
situation in the theory of nuclear interactions makes it 
difficult to escape the conclusion that the Maxwell- 
Yukawa analogy, however suggestive in many of its 
qualitative predictions, is inadequate and that in the 
region of small distances (presumably starting at ranges 
of the order of the nucleon Compton wave-length) 
novel theoretical features must be anticipated. 

It has often been suggested that the typical divergent 
features of present theories may well be due to over¬ 
looking intimate relations between elementary particles 
of various kinds. 6 However this may be, the apparent 
lack of connections between such particles as well as 
between the various dimensionless constants repre¬ 
sentative of interaction strengths calls for more compact 
methods of formulation, with the ultimate ideal of the 
predictability of particles from deeper lying principles. 

Now, whether or not the answer to such questions 
will eventually lead us outside the domain of concepts 
embodied by relativity and complementarity, it seems 
worth while at this stage to ask for models which lead 
to convergent answers. An apparently appealing model 
of this type can be obtained as follows. 

The electromagnetic as well as the meson field 
equations are generally of the prototype 

(□ -K?)'P=P, ( 1 ) 

where p is the source creating a field described by 
Depending on the transformation properties of $ and p 
one' has scalar, vector fields, etc.; k = 0 corresponds to 
photon, and kt^O to meson fields. Finally, the reality 
properties of ^ determine the charged or neutral char¬ 
acter of the field. Now consider instead of (1) an 
equation of the prototype 

F(ni)4'=p ( 2 ) 


6 See attempts to eliminate divergences by compensation: A. 
Pais, Phys. Rev. 68, 227 (1945); Verh. Kon. Ac. Wetenschappen, 
Amsterdam 19, 1 (1947); S. Sakata and O. Hara, Prog. Theor. 
Phys. 2, 30, 145 (1947). 
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*0=na-K< 2 ). (3) 

t=l 

That this is a model of the type we want is clear. First, 
(2) describes quanta of rest mass 7 Ki and thus in a 
trivial way unites particles of various kinds. Second, 
these particles are all bound with the same strength to 
the source p. Third, an equation like (2) leads to a less 
divergent behavior of \p at points near the source p 
than (1) does: take, for example, for p a point source 
with strength g: p=g8(x). Then the static potential 
following from (2) and defined by 


is 


F(A)t=g8(x) 


'Kr) = 



sin kr 


kdk -, 

n-p) 


where r is the distance from the source. The strong 
increase of F with increasing k leads to a decrease of \p 
for small r . In fact, taking in (3) N=2, /cit^k 2 , we have 

^W= (?/4irr)[exp(-)dr)-exp(-K8r)] (4) 

which is already singularity free at the origin. 

It is the aim of this paper to develop the quantum 
field theory of equations of the type (2). Clearly Eq. 
(2) can, in the absence of sources, be derived from the 
following action integral 


L= J 


( 5 ) 


the integration being extended over the whole four¬ 
dimensional space-time (volume element d&). Equiva¬ 
lently (5) can be written as 





£(x„) = J’ V(*) e(x—x / )yp(x')d 4 x / 


where 


1 r 

efe) = -——J dik exp 


( 6 ) 

(7) 


is essentially the Fourier transform of F. For a general 
F containing arbitrarily high powers of Q—we will in 
fact admit that N in (3) may tend to infinity provided 
the infinite product thus arising is mathematically well 
defined—it is readily seen from (6) that the Lagrangian 
density £(xn) will generally depend on the field variables 
at finite distances xj—x^ from the world point under 
consideration. Following a terminology introduced by 
Dirac 8 we can thus say that the present type of theory 
is characterized by a non-localized action. 9 The rela¬ 


7 Throughout this work we put h—c— 1. 

»P. A. M. Dirac, Phys. Rev. 73, 1092 (1948). 

9 The function \(/{x) itself is, however, still perfectly localized. 

The recent suggestion of Yukawa [Phys. Rev. 76, 300 (1949); 


tivistic invariance of L guarantees the existence of a 
symmetric divergence free energy momentum tensor. 10 
Thus, energy and momentum densities can be con¬ 
structed which, like £, have non-localized character¬ 
istics. 

Various investigations have been made of theories of 
this kind from a classical point of view. In particular, 
the work of Bopp 11 should be noted; this author studied 
especially the following generalization of the electro¬ 
magnetic theory: 12 

~j, (8) 

and showed that it leads to a classical theory which is 
not only singularity free and goes over into the con¬ 
ventional theory for distances »k -1 (k can be chosen 
conveniently large), but also does not exhibit the 
so-called runaway solutions which so often occur in 
attempts to eliminate classical infinities. 13 Thus if one 
considers, for instance, a harmonically bound particle 
to be the source of the field the net radiation 
reaction on the particle always leads to a damping, 
notwithstanding the fact that the radiation field corre¬ 
sponding to (8) consists of an electromagnetic field and 
of a mesonic field, the latter having the disagreeable 
property of being negative definite. One would expect 
that such negative energy radiation would somehow 
lead to trouble and it has indeed been shown by 
Feynman 14 that situations may arise in which this is 
the case. 

The quantum-mechanical radiation field (8) consists 
of an assembly of photons of positive definite and of 
neutral vector mesons with a negative definite energy. 
The negative quanta lead to even graver difficulties in 
quantum theory than in classical theory. 15 In Section 
III-A-2 we will discuss this aspect for the general Eqs. 
(2) and (3), and it will turn out that while formally 
any such equation leads to a finite self-energy for the 
electron in the hole theory as long as N> 1, the non-defi¬ 
niteness makes such equations inacceptable. The same 


77, 219 (1950)] to consider non-localized fields, i.e., systems in 
which \p(x) itself is replaced by a non-localized entity, thus 
constitutes a departure of a different nature from equations like 
(1); see however M. Fierz, Phys. Rev. 78, 183 (1950). 

10 See F. Bopp, Zeits. f. Naturforschung 1, 237 (1946); T. 
Chang, Proc. Camb. Phil. Soc. 42, 132 (1946); 44, 76 (1948); 
J. deWet, Proc. Camb. Phil. Soc. 44, 546 (1948); Proc. Roy. Soc. 
A195, 365 (1949); H. S. Green, Proc. Roy. Soc. A197, 73 (1949). 

11 F. Bopp, Ann. d. Physik 38, 345 (1940); 42, 572 (1943); 
Zeits. f. Naturforschung 1, 53 (1946). 

12 This equation has also been proposed by A. Land6 and L. 
Thomas, Phys. Rev. 60, 121, 514 (1940); 65, 175 (1944); by B. 
Podolski et ah, Phys. Rev. 62, 68 (1942); 65, 228 (1944); by D. 
Montogomery, Phys. Rev. 69, 117 (1947); and by A. Green, 
Phys. Rev. 72, 628 (1947). See also B. Podolski and P. Schwed, 
Rev. Mod. Phys. 20, 40 (1948). Recently, a discussion along 
similar lines has been given by L. deBroglie, Comptes Rendus 
229, 157, 269, 401 (1949). 

13 See P. A. M. Dirac, Proc. Roy. Soc. A167, 148 (1938); W. 
VVessel, Ann. d. Physik 43, 565 (1943). 

14 R. P. Feynman, Phys. Rev. 76, 939 (1948); see especially 
p. 945. 

16 See A. Pais, reference 6, Chapter II. §7; P. T. Matthews, 
Proc. Camb. Phil. Soc. 45, 441 (1949). 
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holds mutatis mutandis for “multi-meson equations” of 
the types (2) and (3) involving coupling with nucleons 
(see III-A-2). 

The final results for self-energies on the basis of (2) 
and (3) are essentially equivalent to those obtained by 
Pauli and Villars 16 by means of the “regulator-pro¬ 
cedure” with the single difference however that the 
“regulator condition” for the self-energy here follows 
automatically (see III-A-2). 

It should also be noted that (see III-A-4) the negative 
energy difficulties remain the same whether all in (3) 
are different or whether some of them are equal. 
Thus, for example, the equation 

(□ (9) 

recently discussed by Bhabha, 17 which leads to an 
exponential instead of a Yukawa static potential, still 
has the non-definiteness as a stumbling block. Nor will 
it turn out to lead to any improvement in this respect 
to take some k; to be complex. In fact additional 
complications now occur as will be explained in Section 
III-A-3. 

It may be pointed out here that the occurrence of 
such energies can only constitute a difficulty in the case 
of fields of the Bose-Einstein type (of any integer spin). 
Generalizations of the Dirac equation to equations like 

d^d/dxp ( 10 ) 

y-i 

do not lead to such complications (III-A-5) because 
here an appropriate definition of the vacuum, taking 
into account the exclusion principle, can help out. We 
have investigated whether the use of (10) can be of 
help in eliminating the divergences occurring in the 
charge renormalization of the electron. Such an attempt 
was suggested by the “regularization” of the charge 
renormalization by means of auxiliary fields of the 
Fermi-Dirac type. 18 It will be shown, however, that 
the use of (10) sheds no new light on this problem. It 
seems, in fact, that the construction of a model which 
might conceivably correspond to reality and which 
leads to a finite vacuum polarization is much more 
difficult than finding a model yielding finite self¬ 
energies (cf. also III-B-4). 

As a next step in the discussion of “model theories,” 
one might try to study equations like (2) where F is a 
general integral function of Q On the basis of Weier- 
strass’ product theorem, this means that we now have 
to envisage the presence of exponential functions of the 
dalembertian. It will be shown in Section III-B-1 that 
the new features we now have to consider are all 
exhibited by an equation of the type 

e /(a) ([ |— k 2 )\1/= p. (11) 

16 W. Pauli and F. Villars, Rev. Mod. Phys. 21, 434 (1949). 

17 H. Bhabha, Phys. Rev. 77, 665 (1950). We are indebted to 
Professor Bhabha for communicating his results. 

18 See reference 16 and also D. Feldman, reference 3. 


A special case of this, with /Q)= —[ ], has been 
proposed and investigated by Born and Green within 
the framework of reciprocity theory. 19 Here one encoun¬ 
ters problems of a new type; viz., the lack of “propaga¬ 
tion character.” By this the following is meant: consider 
a space-time point P and its light cone. Now, on the 
basis of our present conceptions, physically meaningful 
differential and integral equations must be such that 
knowledge of any physical quantity in P must be 
obtainable from these equations by specifying sufficient 
initial conditions inside the past light cone only. It will 
be shown in Section III-B-2 that this condition is 
satisfied if F (□) in (2) is a polynomial in □ with 
arbitrary constant coefficients. On the other hand, we 
shall see in III-B-3 that (11) does not have propagation 
character whenever /O)^0. However, the situation 
is radically different for the cases that / is an odd or is 
an even function of its argument. In the latter case 
deviations from orthodox theory show the suggestive 
feature that they seem to average out over space-time 
regions of the order of a universal length X, the magni¬ 
tude of which so far remains arbitrary. 19a Some aspects 
of the situation in this case are discussed in III-B-4 
and in the concluding Section IV. 

As in conventional problems, we shall first study 
uncoupled systems. Now, in the same sense as the 
equation (Q— K 2 )\p=0 has the harmonic oscillator as 
its mechanical model, knowledge of the properties of 
the latter essentially determining all characteristics of 
the former, so the general equation FQ)^= 0 will have 
its mechanical model too. We shall discuss these models 
in the next section. 

II. THE MECHANICAL MODELS 
A. Equations of Motion of Finite Order 

1. Introductory Remarks 

We consider one-dimensional mechanical systems 
(coordinate q ) whose equation of motion, although 
involving time derivatives higher than the second, is 
still of finite order 2N so that it can be written as 

F(D)q= 0, D=d/dt (12) 

where F is a polynomial of degree 2 N. We shall deal 
exclusively with reversible motions so that F is an 
even function of its argument. In this case (12) -is 
derivable from an action principle with the Lagrangian 

L= —qF(D)q. (13) 


19 For a survey and references see M. Born, Rev. Mod. Phys. 
21, 463 (1949). We have been greatly stimulated by this article 
which we had an opportunity to see prior to publication. Also, 
we are much indebted to Dr. H. S. Green for many useful dis¬ 
cussions on this subject. 

19a In this respect the theory shows relationships with the 
modifications of classical electrodynamics as proposed by R. 
Peierls and H. McManus, Proc. Roy. Soc. A195, 323 (1948); see 
also J. Irving, Proc. Phys. Soc. London A62, 780 (1949). 
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In conventional mechanical problems the transition 
to the quantum theory is made by first putting the 
classical equations into Hamiltonian form. We shall 
show in what follows that the quantization by means 
of Hamiltonian methods is always possible in case F in 
Eq. (12) is of finite degree. The first question thus is to 
hamiltonize the system described by (13) and we must 
ask especially how this can be done in a way most 
suitable for subsequent quantization. 

A procedure for deriving a Hamiltonian correspond¬ 
ing to (13) was given long ago by Ostrogradski 20 in case 
F is a polynomial: One defines quantities Qi, Pi by the 
relations 


Qi = D^q, ). 

— 1, 

P t -=5Z,/5(Z?fy),l 


*• ,N , 


(14) 


where D N q is the highest derivative of q occurring in 
L ; 21 8L/8x denotes a variational derivative: 


8L dL dL dL 

—=- D -f-Z> 2 - 

8x dx d(Dx) d(D 2 x) 


Then the Hamiltonian is 

H=P 1 Q i +P&+ • • •+ Pn-iQn+ Pn (DQn) — L } (15) 

where DQn is to be expressed in terms of the P’s and 
Q’ s by means of 21 Eq. (14) for Pn . The Hamiltonian 
equations of motion derived from (15) reduce to (12) 
and a number of identities. 

It is now possible to quantize (15) in the usual way, 
but this method meets with two objections : 

(1) Since the method works only for finite N one can 
at best hope to master the quantization problem for 
transcendental F by approximating F by a finite number 
of terms of a power series, the development of F around 
the origin. In the methods to be discussed approxima¬ 
tions of this kind can actually always be avoided. 

(2) Even for finite N the method is clumsy as can be 
seen, for example, from the non-separable character of 
the Schrodinger equation in N variables corresponding 
to (15). 

One expects however that for finite polynomials the 
problem should be separable since by writing 22 

F=mi+D*/uV, (16) 

i=l 


20 M. Ostrogradski, Memoires sur les equations differentielles 
relatives au probleme des isopirimetres , Mem. Ac. St. Petersbourg, 
VI 4, 385 (1850). See also Whittaker, Analytical Dynamics (1937), 
(Cambridge University Press, London, 1937), fourth edition, 
p. 265. 

21 When using Ostrogradski’s method, it is most convenient to 
perform partial time differentiations on L such that it contains 
D N q as highest derivative instead of D m q, as is the case for (13). 
In this alternative form for L, D N q will occur quadratically which 
ensures that DQn can be eliminated from the Hamiltonian (15), 
while also there will then be no constraints between coordinates 
and momenta, in contradistinction to what would happen if one 
were to start from L in its form (13). 

22 This is always possible apart from a trivial multiplicative 
constant. We will never have to deal with any being zero. 


the solution of (12) is a linear combination of oscillations 
with frequencies At least for the case that all co; are 
real and distinct one would expect the Hamiltonian to 
be some linear combination of oscillator Hamiltonians 
which would make the quantization trivial. Actually 
we will see that the restriction on the co; just mentioned 
is in no way necessary to achieve separability. Thus 
instead of following the Ostrogradski method, which is 
essentially based on a power series representation of F, 
we shall use only product representations and then 
have only to distinguish between the cases of real or 
complex, single or multiple frequencies w*. 

2. Real and Distinct Frequencies 

We start from the Lagrangian (13) where F is given 
by (16) and the co* are real and distinct. It is natural 
to define N coordinates Qi by 

&=n'(i+£ 2 A>/)? (i 7 ) 

/=i 

where the prime on the product sign denotes that the 
ith factor is missing. According to (12) and (16) the 
Qi satisfy the equation 

(D 2 +a>*)Qi= 0, (18) 

while there are clearly no constraints among the Qi as 
their number just corresponds to the number of 
independent solutions of (12). 

Now instead of describing our system by the La¬ 
grangian (13) we try to use another form: 

(19) 

3—1 

We are entitled to do so if we can fix the N constants 
rjj in such a way that L differs from L by at most a 
time-derivative of a function of q and its derivatives. 
We shall show that this is possible, and that the rjj are 
uniquely determined. 

Denoting by the symbol ~ an equality apart from 
an additive time-derivative, and using (16) and (17), 
we have in fact: 


k 3 


VkWi r 

qF*(D)Z - —q, 

k 1 + D 2 /uk 2 


which is equal to L provided that 


V kWkT 


k l+D 2 /a> k 2 F(D) 


( 20 ) 


As the roots of F(D ) are simple we can perform a 
partial fraction decomposition of F~ l with constant 
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(D-independent) numerators. Hence it follows that 

n*=l /«* 4 ^(-«**), 

where (21) 

F'(~u k *) = (dF/d(D>)) D ^ k \ 

Thus we have found the t] k and thereby have obtained 
a Lagrangian L dynamically equivalent to L. 

It is important to note that as F is a single-valued 
function, the quantities F'(—o) k 2 ), and therefore r] k , 
alternate in sign. Thus if we normalize F so that 
F'(—g) i 2 ) is positive, then all r) k (k odd) are positive 
and all rj k (k even) are negative. 

For future reference we mention two further types of 
relations involving the rj k . In the first place one verifies 
that q is connected with the Qj by 

iWQi- (22) 

J-l 

Second, there exist N simple identities between the y} } 
which are easily obtained from (20): 

J2ykW k 2 =l 7 (23) 

i 

ZW n =0, n=2, •••,A T . (24) 

i 

The first of these follows from (20) by putting D— 0, 
while the others are proved by developing both sides 
of (20) in a series in D~ 2 . Next, putting L in the form 

we can immediately write down the Hamiltonian 

ff=LC(iV/4%)+WQ/]- (25) 

i 

Performing the contact transformation 

P r ->P i -(2|, J -|)», Qr^Qri 2|uy|)-* (26) 

and remembering the sign properties of the rjj mentioned 
above, Eq. (25) becomes, in these new variables, 

( 27 ) 

;=i 

This shows that our system consists of a linear combi¬ 
nation of harmonic oscillators, some with positive- 
definite energy (j odd) and some with negative definite 
energy (j even). Quantization yields the eigenvalues 

%=0,1, .... (28) 

l 

Clearly, the “oscillator coordinates” Qi given by (17) 
and their conjugate momenta ought to be related to 
the “Ostrogradski variables” (3) for the corresponding 


problem by a linear contact transformation. This can 
be shown in the usual way by verifying that the 
difference J^PidQi for the oscillator and for the Ostro¬ 
gradski variables is a perfect differential. 

Finally, it may be noted that whereas we can write 
down an L as well as an H in terms of the oscillator 
variables, we have only a Hamiltonian, Eq. (15), in 
terms of the Ostrogradski variables. An Ostrogradski 
Lagrangian does not exist, since (15) is linear in the 
momenta. 

3. Real and/or Complex Distinct Frequencies 

If we admit some of the oh in (16) to be complex, 
the reality of quantities like L and H necessitates the 
occurrence of complex conjugate pairs of frequencies. 
Irrespective of reality properties the formal introduction 
of the Qi by means of (17) is still possible with the Qi 
satisfying (18). Likewise, one can define an L by (19) 
with r) k given by (21). As the problem is still separable 
in the various Qi, ru will be real if oh is real, while if oj m 
and co n are complex conjugates the same is true for rj m 
and rj n . The Hamiltonian is again given by (25). We 
will now pick from (25) a pair of terms corresponding 
to a pair of “complex conjugate oscillators.” Clearly 
then, if we know how to quantize the Hamiltonian 
describing such a pair we will have solved the question 
of quantizing the entire Hamiltonian (25). 

Let j= 1, 2 in (25) denote such a pair. We thus 
consider the Hamiltonian 

H= L(Pmrn)+rn^Qn+ [(iV/4772)+ 

o) = on = o)2*, rj2=rji*, (29) 

We consider first some classical features of our system. 
From (29) it follows that 

(D 2 ~Foo 2 )Qi— (D 2 +o>* 2 )Q 2 =0 (29a) 

as expected. Further, the part of q corresponding to Qi 
and Q 2 is according to (22) given by 

q=Viu 2 Qi+V 2 U* 2 Q 2 ; (D 2 +o> 2 )(D 2 +o>* 2 )q=0. (29b) 

Consider a solution of (29b) of the form 

q—Rie~ at cos(vt+ 0 i)-\-R‘ 2 e at cos(j vt+ 82 ), (29c) 

where we have put 

tu= v-\-ia. 

In contrast to the case of two real roots the motion 
(29c) is in general unbounded in time. Moreover, 
considering the two terms on the right side of (29c) as 
two rectangular coordinates in a plane, the orbit will 
not densely cover a portion of the plane, as in the case 
of the Lissajous motion with incommensurable real 
frequencies, and therefore the motion is degenerate in 
the classical sense. 

With the help of (29a, b) we can express Qi and Q 2 in 
terms of the constants jRi, R 2 , 6 1 , 0 2 and then calculate 
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the energy using (29). The result is (we put 7?i= rj exp i(p) 

H~( 2 RiR 2 /r})[(v 2 —a 2 ) cos(0 1 —0 2 — <p) 

— 2 pasin( 0 i— 02 —<p)J. (30) 

Again in contrast to the real frequencies the energy 
depends on the phases and can thus have either sign. 
In the quantum theory we will find the counterparts 
of the unboundedness of the motion and the indefinite¬ 
ness of the energy. 

Performing the contact transformation 2,3 
Pi — >P 1 (2t]i)\ P 2 -*P 2 (2rn*)* 

Qi—*Qi/(2vi)K Q*->Qt/(2 Vi *)K 

Eq. (29) becomes 


Here n can be any positive or negative integer or zero, 
while X is continuous and ranges from — oo to + . 

The energy spectrum is therefore indefinite, continuous, 
and each level is infinitely degenerate. „,x has been 
normalized in the X-scale. 25 

Let us now consider the matrix elements for the 
operators Q i, Q 2j Pi, P 2 occurring in (31). We have, 
for example, 

SO 

(m, \ f \Qi\n, X") =- [— 1) 

47T(A 


H^P^+JQ i 2 )+|(P 2 2 +co* 2 02 2 ). (31) 

For purposes of quantization it is important to note 
that here the P’s and Q ’s are non-Hermitian operators. 
In fact, taking q in (29b) to be Hermitian, we must have 

Qf-Q*, P 1 W 2 . (31a) 

One can go over to Hermitian variables by the canonical 
transformation 24 

Pi= 

P 2 = W h lp ~ iq ) ~ i(P+ iQ) ], 

i (32) 

Qi= --[ (P~iq)+i(P+iQ)l 

2iu* 

Q 2 = —~[ (p -+ iq) — i{P— iQ )]. 


In terms of these variables H becomes 


H=- oc(PQ+pq) - v(Pq - pQ ), 

where the first bracked has to be properly symmetrized 
to make H Hermitian. Thus the Schrodinger equation is 


- a ( Q -b q— ) + V ( Q - q — ^ 1 = ( iE + a ) V . 

\ dQ dqf V dq dQ/\ 


Putting 

q—r cos^, Q = r sin0, 

we get 

otr{d / dr) — v(d / d0)^p& — (iE-\~a)ty 
with eigenvalues 


E n> \ = nv-\-\a 


and wave functions 


(33a) 


*»,x= {\/2ir)e~ ine ’ (l/r)e ainr . (33b) 


23 This is not a real transformation, but that does not matter 
as the dynamical variables in (29) are not real to begin with. 

24 This transformation is essentially a transition to Fock-like 
variables [see P. Dirac, Quantum Mechanics (Oxford University 
Press, New York, 1949), third edition, Oxford, 1949, p. 136J 
with a subsequent splitting in Hermitian and anti-Hermitian 
parts. 


X 


r 

I dzexp[i(\'—A)z—z]-{-d(M — n—l) 

~ —QQ 

X J dz expp(X'— X)z+s]J, (34) 


where 2 = lnr. The z integrals are badly divergent. Thus, 
corresponding to the classically observed unbounded 
motion the representation for Qi (and likewise for the 
other variables) is unbounded. This does not mean that 
the matrix elements of all functions of p and q are 
unbounded, as is already clear from the finite expression 
(33) for the energy. Indeed, one verifies also that 

(w, V|i(Pi 2 +«^i 2 )|», X") 

= %co(n— i\ rr )h{m— n)8(A r — X"), 

(«, X , |i(P 2 2 +«*^ 2 2 )|», X") 

= |a)*(n+iX // )5(w-w)3(X'-X /, ) > 

so that the energy of each of the complex oscillators 
separately is also finite, although complex. 

In conclusion, it may be pointed out that for the 
reality of H the relations (31a) are not strictly neces¬ 
sary ; one could also have taken 

Qj=-Q 2 , Pf=-Pt, (31b) 

corresponding to anti-Hermitian q. In the case of a real 
oscillator, one chooses either a Hermitian or an anti- 
Hermitian q in order to obtain a (positive or negative) 
definite energy; the same holds for the case of the 
distinct real frequencies discussed in the previous 
section. For complex frequencies, however, each choice 
leads to an energy which is indefinite to begin with, and 
therefore (31a) and (31b) cannot be distinguished from 
each other. 


4. Presence of Multiple Frequencies 

We first consider the case that only a double real 
frequency is present, so that L is 

_ L=-q(D 2 +m 2 ) 2 q. 

26 This is meant in the customary sense of the continuous 
spectrum, i.e., one takes m=n and integrates over X" from X'—e 
to X'+e where e is a small quantity. 
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The general solution of the equation of motion is 

q = Ri cos(a)oH-0i)+-^ cos(a)o^ - f* 02 ) J (-35) 

which in general is again not bounded and degenerate. 
For the energy one finds 26 

4cuo 2 2^2D^iwo sin(0 i— 02 ) — ^ 2 ]] 

which as in the complex case is indefinite. For the 
Hamiltonization process the partial fraction method of 
the previous two sections is clearly inapplicable, while 
the Ostrogradski method though feasible is again incon¬ 
venient. It is natural, however, to try here to consider 
the double frequency as the confluence of two distinct 
real frequencies. For the latter one has 

L=-?(D 2 +u 2 )(D 2 W)<7 
«[Q^+c^Ci-eoP 2 w)eo]/(« 2 -w), (36) 

with 

Q,= {D^+^)q, Qo=0?+'f)q, 
q= (0i—Qo)/(a> 2 —o> 0 2 ). 

The occurrence of a factor (w 2 —a> 0 2 ) -1 in (36) indicates 
already that the transition w—>w 0 has to be performed 
with some care. We put 

oj = coo+e 

where e is considered infinitesimally small, and also 
Qq—Qv\- *(?2. 

Then L becomes in the limit e—>0 

« 0! 2 + (1/COO) {DQl • DQt-^QlQt), (37) 
and in the s?me limit 

q= —Q 2 /2oi Q . 

One easily verifies that (37) leads to the same general 
motion (35) for q, which justifies the procedure. From 
(37) one finds for the Hamiltonian 


with n a (positive, negative, or zero) integer while r f 
ranges continuously from 0 to 00 . Thus, as in the 
complex case, the energy spectrum is continuous, 
indefinite and each level is infinitely degenerate. 

The general procedure for dealing with a mixture of 
single and multiple frequencies of arbitrary multiplicity 
is now clear: dissolve first each w-fold root co into n 
neighboring ones 03+e i} i= 1, •••,», with e^e;. Then 
develop in L with respect to the e* and pass to the limit. 
As for our further discussions a detailed analysis of 
this procedure seems hardly necessary, we will not 
pursue this question any further. However some addi¬ 
tional comments regarding multiple frequencies will be 
found in the Appendix. 

To recapitulate the results of Section A: if L contains 
more than one frequency the energy is always indefinite, 
both classically and quantum mechanically. For the 
case of real frequencies we saw that this follows immedi¬ 
ately from the alternating signs of the rj k . In the other 
cases this simple argument cannot be taken over, but 
from the analysis of the Schrodinger problems it was 
found that the same conclusion holds. It is this indefi¬ 
niteness which will constitute a grave difficulty in field 
theoretical applications. 

B. Equations of Motion of Infinite Order 

1. Introductory Remarks 

If F(D) contains derivatives of infinite order we have 
essentially to do with integral equations. It is therefore 
often more convenient in this case to use instead of 
(12) the following form for the equations of motion 

fK(t-t')q(t'W= 0, (38) 

where the time reversibility is now expressed by the 
requirement of the evenness of K as function of its 
argument. Corresponding to (38) L has the form 


H —coo (P1P2+ Q1Q2) — < 2 i 2 

which, on performing the canonical transformation, 


.= -?(<)/*(/■ 


t')q(t')dt'. 


(39) 


Qi~*Qh Pr-^Pi—Qi/^oiQj 

Qi~>— iVF(?i/2a>o, P 2 —>P 2 , 

becomes 

^ ==CJ o(P 1 0 2 -P 2 Ql)-KQl 2 +G2 2 ). 

Put 

Qi = r cos0, Q 2 =r sin0, 
then the Schrodinger equation is 

d*/d0=(£+|r 2 )*, 

so that 

E=no>-r f 2 /2, 0) = 8(r-/)e in9 /(2irr')K 


The connection between (12) and (13) on the one hand 
and (38) and (39) on the other is expressed by 


K(i) = J dke ik ‘F(ik). 


(40) 


We will restrict ourselves to functions F which are 
integral functions in the sense of complex function 
theory. It is known that according to Weierstrass such 
functions can be represented in the form 

F(D) = e wmi+D 2 M, (41) 

i 


26 As can, e.g., be shown from v he Ostrogradski form of the 
Hamiltonian. 


provided converges. Equation (41) expresses F 

in terms of its zeros and of an exponential factor;/(D) 
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is again an (even) integral function determining the 
behavior of F at the point at infinity. The motivation 
for the choice of functions of the type (41) is that in 
corresponding field problems the zeros will each corre¬ 
spond to quanta of a certain mass whereas the expo¬ 
nential is suggestive of a cut-off factor. 

Again the question arises of the quantization. It is 
clear that if no exponential is present (i.e., /= 0), the 
methods of the previous chapter can be applied immedi¬ 
ately. We illustrate this with a simple example: If 
F(D) = coshaZ) we have the product representation 

FM-mi+PM), OH=Tc(2i— l)/2a. 

l 

Thus in our terminology we have an infinite number of 
distinct real frequencies. Now we can define an infinite 
number of Q z by using (17) (for N= qo) and an L as in 
(19). From (20) one sees that the rj k are to be determined 
from the relation 

E 7 ?* &k 2 / (l-\-D 2 /o) k 2 ) = l/coshaD. 


tonian becomes awkward since one would always have 
to approximate the exponential either by a power series 
of finite order (in the Ostrogradski method) or by a 
finite product representation (in the oscillator method), 
and then in some way have to pass to the limit. How¬ 
ever, one can circumvent the use of the Hamiltonian 
for quantization by following a method advocated 
especially by Heisenberg. 29 The idea here is to start 
from (1) the equations of motion, (2) the expression of 
the energy E in terms of q and its time derivatives, 
(3) the quantum-mechanical definition of time deriva¬ 
tion which follows from attributing to E the role of 
time-displacement operator, thus 

[E., v {q,q,q, ••■)]= ~i<P (43) 

for any function ip. From these three requirements the 
eigenvalues of the energy must be determined and the 
representation of q must be found. 

We will now apply Heisenberg’s method to the 
system (42). First of all, we need an expression for the 
energy, to the derivation of which we now turn. 


The partial-fraction series for cosh r l aD is well known 
to be 27 


so 


4 * (— l) i_1 
sech aD=-E 


Vk- 


X 1 2k— 1 l + £ 2 /W 2 
4 (— l) fc_1 16a 2 (-1)*- 1 

x (2k-1W x 3 (2k— l) 3 


2. Classical Theory—General Expression for the Energy 

From the general definition of E as the integral of 
the motion for a time displacement, one readily finds 
that 

E= L D n q-\_bL/b(D n q)~]-L 


One notes again the sign alternation of the rj k . The 
identity (23) is easily verified. 28 

If there is an exponential present in F, the problem 
can still be reduced greatly by the above method of 
decomposition. In fact it is readily verified that for the 
Lagrangian (13) with F given by (41) we can introduce 
an L~L, where 

L=ZrikQke^( l+D 2 /o> k 2 )Q k , 

with the same Q k and rjk as before. Therefore the 
quantization problem will be mastered if we can deal 
with a system whose Lagrangian is 

L= —qFi(D)q, F 1 (Z)) = |^>(D 2 -fco 2 ), (42) 

i.e., for a system with only one frequency and an 
exponential factor. 

For a system of this kind the generalization of the 
various methods described above for finding a Hamil- 

27 See Whittaker-Watson, Modern Analysis (Cambridge Uni¬ 
versity Press, London, 1940), fourth edition, p. 136. Quite gen¬ 
erally the possibility of decomposing F~ l (D ), F given by (41) 
with /== 0, and So of* convergent, in a partial fraction series is a 
consequence of the Mittag-Leffler theorem of complex function 
theory. 

28 The identities corresponding to (24) now are 

00 

2?7£a>& 2n = 0, 2 

l 

which are divergent series but by the known methods can be 
summed to zero. 


= E D n qZ (-\) m D™[_dL/d(D m + n q)~\-L. 

n=l 0 


We will consider the Lagrangian to be of the general 
form (13). Its value then is zero in virtue of the equa¬ 
tions of motion (which are supposed to be given and 
thus may be used freely in contradistinction to the 
situation in the Hamiltonian method). It is convenient 
to introduce the Fourier transforms of the various 
quantities. Writing 

F0D) = I> n Z>", 

o 

q(t) = (l/2r)f R{k)e~ ikt dk, 


m 


= (1/2 ir)J 


S(s)e~ i8t ds , 


one finds 



E E (-1 ) m Xn+ m 

n=l ro=0 


f dk-(-ik) m -{-i(s-k)} n R{k)R(s-k) 

-DO 


29 W. Heisenberg, Zeits. f. Physik 123, 93 (1944). 
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Putting and retaining m, one can carry out 

the summation over m. The subsequent summation 
over p can be expressed in terms of F. The result is 

S(s) = (l/2rs)J dk(k-s)R(k)R(s-k) 

XtF{i(k-s)}-F(ik)l (44) 

and we note that as the equation of motion in 
Fourier form is F(—ik)R(k) = 0 and as F is an even 
function of its argument, s£(s) clearly is zero as a 
consequence of the equation of motion; this expresses 
the energy conservation. For the L given by (42) we 
have 

< 7 — A. cose*?/, R—wA 8(k —co)]. 

One finds that S(s) is proportional to 6(5) and eventually 
that 

E=iAW^\ (45) 

If 

F(D) = i x (D)Il(l+D*M, (46) 

i 


where x is any function without zeros one may con¬ 
veniently write 

*7 =: ^L ]VnA n COSW„/, R— n[]5(^~|“^n)-i“ 8(k C**n)3> 

n n 

where the rjk are the constants given by (20) for 
One then finds 

E=iYiVkx(^k)(*3k 2 A k 2 . (47) 

k 

This additivity of the energy follows immediately from 
the corresponding additivity of the Lagrangian as found 
in the previous section. 


3. Quantum Theory—Diagonalization of the Energy 

We will discuss now the quantization of the system 
described by (42), making use of Fourier transforms 
throughout. Then (43) is clearly satisfied if 

[£(£),£]=££(£). 

In the representation in which E is diagonal this means 
that the matrix elements of R must be of the form 

(£'| R(k) | E") = 8(E'—E"—k)(E'\R\ E"). (48) 

As the equation of motion is Fi(ik)R(k) = 0, it follows 
that R(k) is zero except for so that we may 

write for (48) 

(E' | R(k) | E") 

= 8(k 2 - co 2 ) 8(E' -E"-k) (. E ' \Ro\E”)- (49) 

With this form of R all requirements are fulfilled apart 
from the diagonalization of the energy, the condition 
for which is 


1 

2t 



E")e~ ist =E f b(E f -E tf ) 


where S(s) is given by (44) with F—F 1 . From (44) and 
(49) one readily finds 

1 

FJ =-{| (E f | R 0 1 FJ+ co) | 2 

I67r 2 

+ \(E'\R 0 \E f -u)\ 2 }. (50) 

This equation should determine both the spectrum 
of E and the explicit form of the matrix elements 
(E! I Rq I jE/zbco). In the case of an ordinary oscillator 
(/= 0) a solution clearly is 30 

#=(»'+*)«, 

(Ef | Ro | E") = 2ttv2[(w / co) *«(»' - n"- 1) 

+ (» , '«)*5(» / -n"+l)]. 

For our system with the exponential factor one therefore 
must have 

E'=(n'+$) ue'M 

while (E'\Ro\E") is the same as for the ordinary 
oscillator. Thus q mn) the matrix element of q in the 
energy representation is the same as in the case/=0. 

It should be noted that if / is of the form (Z^+co 2 )*, 
both the eigenvalues of the energy and the matrix 
elements of q are exactly the same as for an ordinary 
oscillator. Hence in this case the exponential factor has 
no effect whatsoever, both classically and quantum 
mechanically, as long as one confines one’s attention to 
the homogeneous problem, i.e., to the case that no 
outside forces act. 

As soon as such forces are present the exponential 
factor has a profound effect, however. Consider, for 
example, the problem of an “exponential oscillator” 
coupled to an ordinary one as described by a Lagrangian 

L= -iq exp[X 2ft ( J D 2 +co 2 ) fc ](Z) 2 +co 2 )g 

— !#o(Z> 2 +o>o 2 )<7o+ e (?o q y 

where X has the dimensions of an inverse frequency. 
The equations of motion are 

exp[X 2fc (Z> 2 +co 2 ) *] • (D 2 -\-co 2 )q= eq 0 , (Z> 2 +aj 0 2 )?o= eq. 

Now solve in the usual way by putting 

q^e ivt , qo^e ivt 

where v must satisfy the condition 

exp[X 2fc (co 2 — v 2 )*] • (co 2 — v 2 ) (coo 2 — v 2 ) = e 2 . 

In contrast to the conventional problem with X = 0 one 
does not find two solutions for v, one near co and one 

30 The uniqueness of the solution for the harmonic oscillator 
has been studied by E. P. Wigner, Phys. Rev. 77, 711 (1950). 
Wigner finds that for this case the results of the present method 
are not necessarily the same as those obtained by the Hamiltonian 
method. It has kindly been pointed out to us by Professor Pauli 
that uniqueness can be obtained by requiring the commutator of 
q and $ to be independent of the particular choice of E (a result 
already implicit in the last paragraph of Wigner’s paper). We 
are indebted to Professors Wigner and Pauli for interesting 
discussions on this topic. 
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near cjo, but rather an infinity of solutions, 31 so that as 
soon as € 5^0 the motion becomes a great deal more 
complicated. In field theoretical applications this dis¬ 
tinction between the homogeneous and the inhomo¬ 
geneous problem will turn out to be of great importance. 

in. FIELD THEORETICAL APPLICATIONS 


introducing the following set of field variables 

Aw=f<n)A» 

A^=n'(i-n/rfHn/», s )A„ *= 1 . 2 , - - - 

With the convention p 0 =0 and using partial space-time 
differentiations we can replace (52) by 


A. Absence of Exponentials 

1. Introductory Remarks 


£«- 



dA^ 
dx v dx v 


\i?A^ 


) 


With the help of the methods developed in Section II 
we are now in a position to discuss a number of questions 
concerning field equations of the general type (2). For 
the function F we first limit our choice to expressions 
like (3) (with a finite or an infinite number of factors) 
and defer until Section III-B the discussion of more 
general types of integral functions F which include 
exponential functions of the dalembertian [3 
We begin with the case that \p in (2) is the electro¬ 
magnetic field and the p is the current vector due to the 
presence of electrons, thus 

/0)04,=-i,, /o=iki-qw), (si) 

i 

where jp is the conventional Dirac current operator. 
It is thereafter pretty straightforward to see what will 
happen for other choices of \p and p in (2). Some 
comments on the application to the coupling of mesons 
with nucleons are made in the following subsection. 

In an analogous manner to the developments in 
Section II we shall make increasingly general assump¬ 
tions about the p*: at first that they are real and 
distinct and later that there are also pairs of complex 
conjugate p; and multiple occurring p t *. In subsection 5 
we discuss generalized equations of the Dirac-type (10). 
We now start with the discussion of (51) for real m. 


2 . Real Masses—The Self-Energy of the Electron 

The Lagrangian density from which the generalized 
electromagnetic Eqs. (51) as well as the wave equation 
for the electron can be derived is 


1 dA^ dAp 

£ =-/(□)- jnA„-$ (tA+ 

2 dx v dx v 


(52) 


j»=**&/nf- 


(54a) 

k 

with 

fo=l, f*=l/p* 2 /W); *= 1 , 2 , •••. (54b) 

The transformation from (52) to (54a, b) has been made 
by following step by step the derivation of (19) from 
(13). The coupling term in (54a) is obtained from that 
in (52) by using the analog of (22); i.e., 

As to the derivation of the values (54b) for the f*, one 
readily verifies that the determining equation (which 
here plays the same role as (20) for the mechanical 
model) now is 

i/ana --ow)=ovo -ejv a - □/«*). 

i % 

In particular one finds f 0 = 1 by multiplying by □ and 
then putting □=0. The identities which are the 
counterpart of (24) now turn out to be 

N 

0, n=0, 1, • • •, 2V—T. (55) 

i=0 

Also the sign alternation of the constants is of 
course found here again as in II-A-2. We have normal¬ 
ized our definitions of the A^ k) so that f 0 =l which 
means that the part of (54a) which contains ^4 M (0) has 
precisely the structure of the standard electromagnetic 
field Lagrangian with interaction, so that we now 
identify A^® with the electromagnetic potential. 
Similarly, for k^O, A p w is now a neutral vector meson 
field with meson mass p^. Following the reasoning 
which led up to (28) one sees that for even (odd) k the 
field energy is positive (negative) definite. 

From (54) it follows that, for all k , 


Here yp is the wave function of the electron, yp its 
adjoint; k is the electron mass, e its charge. We impose 
on A\ the supplementary condition 

dA /i /dx fJi =0. (53) 


which means that each field is coupled to the electron 
with the same strength e. Further it follows from (53) 
and from the independence of the fields that 32 


Following the procedure of II-A-2, <£ can be split by 


dAjn/dx^ 0, all A (56) 


31 This is an immediate consequence of the Picard theorem. 
Roughly speaking, one might say that for e = 0 an infinity of roots 
have collapsed to the point in infinity but are “drawn back” from 
there as soon as e diners from zero. Clearly, the transition e—*0 
is highly non-uniform. 


32 It is well known that the essential role of relations of the 
type (56) is to insure the definiteness of the energy for the field 
concerned. In view of the fact that (56) is still insufficient for 
obtaining a positive definite energy for the assembly of fields 
under consideration, one may ask whether the supplementary 
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For the quantization of the fields it is important to 
note that the momentum conjugate to A^ k) (x, t) is 
0* Consequently the covariant commutation 
relations are 

A v {m) (x')^=i5p V d km -(l/tk)&k(x—x f ), (57) 

where A k (x) is the well-known four-dimensional “8- 
function” involving the mass p*. 33 To each A* we can 
adjoin functions A* and A* (1) in the manner indicated 
by Schwinger. 33 In addition, we will assume in parallel 
with (57) that the vacuum expectation value of the 
anticommutator of ^4/ fc) (x) and A v {m) (xf) is given by 

= Wfc»‘( l/tk)A k W(x-x'). (58) 

To the interpretation of (58) we will return presently 
and will now first consider the self-energy of the electron 
in the order e 2 as it can be calculated by means of the 
formal relations (57) and (58). 

We note first that it will consist of a sum of terms 
each being the contribution of a separate field A^ (k) . 
Evidently such a term contains a factor $T 2 since we 
deal with an effect of the second order in the interaction 
energy which according to (54a) is weighted with f*. 
Furthermore it follows directly from (57) and (58) that 
the term we consider contains another factor f* -1 . 
Hence the total self-energy will be of the form: 

8m~ H tkW(JJLk ), (59) 

k 

where W(ji k ) is the standard expression for the self¬ 
energy of an electron coupled to a neutral vector meson 
field with meson mass ix k ; therefore 

e 2 K r 1 

8m = — I rf«(l+«)53f w [(C— In?;) 

8ir 2 Jo n 

- In(w 2 * 2 +p n 2 (l -«))], 

where C— 0.577- • •, and tj is a cut-off which ultimately 
tends to zero. Now according to the first of our identities 
(55) we have 

Ern=o (60) 

» 

and we see therefore that this relation guarantees the 
finiteness of 8m . 

Thus theories of this type automatically lead to the 
finiteness of the self-energy but this result has been 
obtained at the cost of an unreasonable assumption. 
In fact the use of (58) for the vacuum expectation value 

condition (53) is the most appropriate one for the problem on 
hand. It would seem most unlikely, however, that the indefinite- 
ness can be avoided in such a way that the ultimate aim of 
obtaining convergent self-energies can still be reached. 

_ 33 The conventions followed in the definition of A*, as well as 
A k and Afc (1) to be introduced below are those of J. Schwinger, 
Phys. Rev. 75, 651 (1949), Eqs. (A.29), (A.21), and (A.37). 

For k^Q one has strictly speaking to add to the right of (57) 
terms involving second derivatives of A. As those terms can 
effectively be transformed away (see P. T. Matthews, Phys. Rev. 
76, 1254 (1949)) we have not bothered writing them down here. 


of functions of the A^ k) implies that for all fields one 
has chosen as vacuum definition 

0, (61) 

where SP is the vacuum state vector and A^ +(k) denotes 
the positive frequency part of A^ k) . For even &, (f *>0), 
this means that no annihilation of quanta of positive 
energy is possible. For odd k , f* is negative so that the 
r61e of creation and annihilation is interchanged. Hence 
for odd k (61) means that no creation of quanta of 
negative energy is possible. This condition formally 
ensures that the vacuum is a state of minimum energy. 
However, a definition of the vacuum which forbids 
creation processes is certainly completely unphysical 
and, since the quanta obey Bose-statistics, a reinterpre¬ 
tation by means of a filling-up of the negative energy 
states is impossible. Alternative interpretations with 
the help of an indefinite metric in Hilbert space are 
equally unsuccessful (see Matthews 15 ). Theories of this 
kind therefore are physically inconsistent. 

The procedure here described for obtaining finite 
self-energies is closely connected with the regularization 
procedures developed by Pauli and Villars 16 and others. 
As has been emphasized by these authors, the use of 
“regulators” is a computational device which amounts 
to making certain divergent expressions finite. An 
essential step in this procedure is the introduction of 
certain parameters, called d and Mi in the cited paper. 
Between these, certain relations must be postulated to 
achieve the desired convergence. On closer comparison 
of the above results with the regulators one readily sees 
that the d and Mi just correspond to our and p*. 
Thus there are two differences between the regulator 
procedure and the models of the present section. 
First, the “auxiliary masses” Mi formally introduced 
in an advanced stage of the calculation, here correspond 
to actual masses p; introduced from the outset and 
linked with fields. Second, the regularization condition 
XX=0 which has to be postulated in reference 16 is 
equivalent to the identity (60) which here follows 
automatically. Thus the formalism in this section may 
be said to be a realization of the regulator type of 
theory, but as we have just seen, it is not an acceptable 
scheme. 

Actually (60) is only one of the many identities which 
are given by (55) (for suitably large N). In the present 
case only the first of the relations (60) is needed for con¬ 
vergence (which amounts to saying that the self-energy 
is finite for any iV> 1), but in related problems further 
similar relations come into play. Thus it is clear that 
the formalism here developed is immediately adaptable 
to the case where ^ in (2) is some meson field and p the 
corresponding nucleon source. Also in this case one will 
have identities closely analogous to (60). Now it is 
well known that for all couplings involving derivatives 
of the meson field potential (like the tensor coupling in 
the vector theory) the self-energy of the nucleon in the 
first non-vanishing order exhibits a quadratic and a 
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logarithmic divergence. We leave it to the reader to 
verify that in such a situation the first two identities of 
the type (60) suffice to obtain convergence. 34 

3. Complex Hi—Unbounded Interactions 

In the previous section we have seen how the occur¬ 
rence of negative energy states for the uncoupled 
system prohibits a logical physical interpretation of a 
formalism based on (52) and (51) with real masses m. 
These negative energy states were encountered already 
in the corresponding real frequency model of II-A-2. 
In II-A-3 it was shown that for complex frequencies 
such negative energies also occur. Now the latter model 
will play the same role for the case of complex u as the 
real frequency model did for real ju. It is therefore 
clear that we will run into the same negative energy 
trouble whether or not we admit complex conjugate 
pairs of Hi to occur in (51). In certain respects the 
situation gets even worse, however, when complex Hi 
are present. It is this qualitative difference between the 
cases of real and of complex m which we want to com¬ 
ment on now. 

Also if there are some pairs of complex conjugate h% 
in (51) we can still write down the Lagrangian density 
(52) and can again perform the decompositions leading 
to (54a, b). Let us first consider the case of the absence 
of interactions (e=0). The field quantization does not 
differ in any essential from that of more conventional 
situations. By performing a spatial Fourier development 
one can show by standard methods that the total 
Hamiltonian referring to a pair of complex conjugated 
M is a linear superposition of pairs of complex conju¬ 
gated oscillators of the type (31): 

^=iECPi 2 (k)+co A 2 0i 2 (k)+P2 2 (k)+^* 2 02 2 (k)], 

k 

where 

co* 2 = M 2 +|k| 2 , (62) 

with complex h- According to (33) the eigenvalues for 
the kth mode are 

nvk-\-\otk, 

where, as in (29c), we put cok— Vk~\riock. According to 
(62) the dispersion law of the modes of the fields is of 
an inherently complex nature so that no particle 
attributes can be ascribed to these modes in the manner 
used for real ju. 

When dealing with the problem of interaction we 
must be prepared for a queer situation. For the inter¬ 
action energy is according to (54a) linear in the A p , 
and hence, as far as the contribution from complex Hi 
is concerned, linear in field amplitudes which according 
to (62) are unbounded as a function of time. It is 
therefore clear that a coupling like yV4 M (c) , where A^ c) 

34 In this case p in (2) is proportional to a coupling constant of 
the dimensions charge X length. This length should be considered 
to be independent of the k x occurring on the left-hand side of (2). 

Added in proof. See also W. Thirring, Phys. Rev. 77, 570 (1950). 


stands for the contribution to A M of the complex Hi, 
cannot be treated as a time dependent perturbation 
notwithstanding the fact that it is proportional to a 
small parameter, because the coupling is not periodic 
in t. If the coupling had only exponentially decreasing 
modes, one might try to use much the same methods 
as are employed in problems of line-breadth and reso¬ 
nance fluorescence. However the modes which increase 
exponentially with time and which are of course the 
real source of complication, make this impossible. 

Rather than discuss in detail the field theoretical 
applications we shall exemplify with the help of a simple 
model what may happen under such circumstances. 

Consider a pair of complex conjugated oscillators 
(31) coupled to an ordinary oscillator in the following 
manner: 

ii=UPi 2 +*Vi 2 )+hW+u*V2 2 ) 

+ 2(^1 o 2 +wo 2 Qo 2 )+g((?i+ Q 2 )Qo, (63) 

where all quantities with the subscript 0 refer to the 
real oscillator; the coupling is proportional to a constant 
g which we may take as small as we like. Hence it 
must be possible to transform away the coupling in H 
by means of a point-transformation in such a way that, 
in terms of the new variables, H represents an uncoupled 
system, again consisting of one real oscillator and a 
pair of complex conjugated oscillators, where the real 
as well as the complex frequencies have undergone a 
small shift from their original values. By expressing P 0 
and Qo in the new variables one can calculate the 
expectation value of the “old” real oscillator energy 
for a given state of the new system. It is then readily 
seen that part of this expectation value comes from the 
diagonal elements of 

const.(P 1 'P 2 '+Q 1 , <2 2 '), (64) 

where the primed quantities are the new variables in 
the sense just explained. Using the eigenfunctions (33a), 
one easily shows that the diagonal elements of (64) are 
infinite. 35 

Infinities of this type have, of course, nothing to do 
with the divergences of present field theories which are 
due to the infinite number of degrees of freedom of the 
fields. They rather remind one of the pathological 
runaway solutions in electron theory which are well 
known to occur if we attribute to the electron an 
infinite sink of mechanical energy. 

Thus the introduction of complex Hi has brought us 
farther from, rather than nearer to, an acceptable field 
theoretical model. We will not pursue this line of attack 
further. 

3. Multiple Real Masses 

From the results of II-A-4 we infer that again the 
negative energy difficulty prevails. We also met in the 

36 All other contributions to the expectation value of the energy 
of the old real oscillator are finite and thus of no interest for the 
argument. 
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corresponding mechanical model with the feature of 
unbounded motion, although the increase of the motion 
with t is not so rapid as in the complex case. The 
multiple masses thus will occupy a position somewhat 
intermediate between the real and the complex case. 
We have not investigated this category in any further 
detail. 

Recapitulating the results of the last three sections 
we may state that any generalization of the electro¬ 
magnetic or mesonic field equations of the kind (2) 
where F is an arbitrary polynomial in Q may lead to 
an elimination of divergences but only at the cost of 
other unnatural features. 

4. Multi-Mass Dirac Equations 

The partial fraction method can also be applied to 
equations of the type (10). Let us start from a 
Lagrangian 

1 pF(— k 0 )F( 7 A) (7 A+ (65) 

with 

F{x)=Y[(\+x/n,). 

7=1 

The constant factor F(—kq) has been introduced for 
reasons of convenience which will presently be clear. 
Define now 

^0 = F(tA)^ 


complications of the previous sections, as we have now 
to do with Fermions. Hence by interchanging the roles 
of particle and antiparticle for the fields with odd k as 
compared with the even ones one has obtained a 
reasonable formulation for an assembly of A 7 +l fields 
of the Fermi-Dirac type. The normalization of (65) 
with a factor F(—k 0 ) led to £ 0 = 1 and we can therefore 
directly identify the term of (66) with j— 0 with the 
Lagrangian for free electrons (mass k 0 ). 

The quantization can now be performed by putting 36 

{¥*(*), ^iMhit-^+V^Atx-x') 

where the index k on S indicates that this quantity 
refers to the mass n k . A factor (—1)* has to be intro¬ 
duced to obtain positive probabilities for such quantities 
as the number of particles in a given volume. The 
main formula for computing vacuum expectation values 
becomes 

^zM])vac= (— t) k+1 £k~ l dkiS k il) ( x — %') (66a) 

where the factor (—1)* expresses the effect of the 
interchange of particles and antiparticles for odd k on 
the definition of the vacuum. 

In introducing the electromagnetic interactbn one 
will first of all try to replace d„ in (65) by 

— dp ^ 

and add the Lagrangian S$ e im of the electromagnetic 
field. We express this by 


*J= ( 1+7 AAo)IT(i+7 i= l, • • •, N, 

k 


and the adjoint quantities accordingly. Now replace 
£ m by 

N 

"If^itTA+Kj)^. (66) 

0 


The & are to be determined from 


at Kj 1-fyAo F(— K 0 ) 

f 0 + £ £/- —■=—• 

1 Kol+y/Kj F(y) 


(67) 


Hence, putting y— — k 0 : 

fo=l, (68a) 

while, after developing F~ l (y) in partial fractions, one 
readily finds that 


kq 2 F(— K 0 ) 

-, 7>0 

Kj 2 (Kj-Ko)F'(—Kj) 


(68b) 


One also verifies that 


[1/ kqF(— ko)]L£*k A. 

From (68a, b) it follows again that £*>0 for even k 
and <0 for odd k. Here of course the occurrence of 
negative does not lead to the negative energy 


£ = £4<V^A*)+£e^. (69) 

Next one may also in the definition of T'* replace d M by 
Dp so that the matter fields would depend implicitly 
on Ap. It is not difficult to see, however that then 

£=£«0 M ->Z> M )+£ez« (70) 

is not equivalent to £ because it leads to different 
electromagnetic field equations than does <£. Therefore 
the decomposition method (which in the case of multi¬ 
mass equations of the Bose-Einstein type also works 
when an interaction is present) here breaks down. This 
arises from the fact that according to (69) the coupling 
would no longer be linear in A It is an open question 
whether a theory based on (69) leads to an improved 
situation with regard to such questions as the vacuum 
polarization. 

Alternatively, one might try to use (70) itself as a 
starting point for the description of a system of spin 
l-fields coupled with the electromagnetic field, i.e., one 
might, so to say, forget the way (66) was obtained 
from (65). There is no objection to this, but no ad¬ 
vantage either. For it is readily seen that in the compu¬ 
tation of the charge renormalization the £* drop out 
entirely; therefore the vacuum polarization is essentially 
unaffected. 

36 We define S(x ) and S (1) (x) in the same way as J. Schwinger, 
see esp. Phys. Rev. 74, 1439 (1948), Eq. (2.29) and Phys. Rev. 
75,651 (1949), Eq. (1.68). 
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B. Field Equations of the Exponential Type 

1. Introductory Remarks 

We now turn to a discussion of the main features of 
equations of the type 

e m) n o - = — P- 

i 

From the remarks made at the beginning of II-B-1, 
it follows that the problem is separable up to the stage 
at which one has to deal with 


g /(D) (| | — k 2 )^— — p . (71) 

In view of applications to electromagnetic phenomena 
we put k= 0 ; none of the essential new features we shall 
meet presently are affected by this specialization. 

Equations of type (71) go so far beyond the partial 
differential equations usually encountered in field 
theory that it seems necessary to investigate first 
whether the basic mathematical problems and results 
pertaining to differential equations of the wave equation 
(hyperbolic) type have their counterpart in the present 
case. It may therefore be well to recapitulate briefly 
what the two fundamental problems and the two 
principal results are for the ordinary wave equation 

04'= ~P- (72) 

(a) The initial value or Cauchy problem. —A solu¬ 
tion of the homogeneous wave equation is sought 
when \p and its normal derivative are given on 
some space-like surface. 

(J3) The source problem .—A solution of the in¬ 
homogeneous equation is asked for, such that p 
and bp/bt are zero at t—~ oo. These initial 
conditions correspond to the requirement that the 
source p exerts a retarded action. 

It is well known that (a) and (j3) have a unique solution 
and that they are intimately related to each other. 
This relation can, in the language of present day field 
theory, be stated in terms of the connection between 
the D- and the ZMunction 37 

D(x) = — i sign (t)-D(x) (73) 

with sign(/) = ±l for 0. The solution of (a) is: 

p(x)= f d<r|~Z)(x—»:')- \p(x') — D(x— x')\ (74) 

J L bn dn J 

where b/bn denotes differentiation normal to the 
surface on which the initial values are specified and 
over which the integral (74) is extended. The solution 
of (JS) is 

1 K*)=fd i x'D Tet (x-x') P (x'). (75) 

37 See e.g., J. Schwinger, Phys. Rev. 75, 651 (1949), Eq. (A.l). 
The explicit representations for D and D given in reference 33, 
Eq. (A.29) (with k 0 =0) and (A.23), are also used here. 


Here D ret is the Green function for the retarded type 
of solution. It is given by 


A 


= —!— fdk fdw 

(2ir) A J J 


expp(kx— 

k 2 —&> 2 


1 

= — 8(t- 
4wr 


r), 


(r= M) 


(76) 


with the specification that in the complex co-plane one 
shall integrate over a path parallel to and above the 
real axis. By integrating below the real axis one gets 
the Green function D a dv for the advanced solution. 
We have 

S(*) = KAe t +£adv); (77) 

all three functions in (77) satisfy 

D0=-«(*). (78) 

After this connection between (a) and (0) a second 
fundamental aspect may be recalled; viz., that for both 
problems the solutions do not depend on all initial 
values or values of the source but only on the values 
in a finite domain of dependence. [In fact for (0) this 
domain is the past light cone, and for ( a) the part of 
the space-like surface bounded by that cone.] This is 
the mathematical expression for the propagation char¬ 
acter of hyperbolic equations. For the wave equation 
(72) the propagation character is physically obvious 
and can also be stated as 

< w > 

i.e.; both D and D vanish outside the light cone. 

The suitable generalization of the problems (a) and 
(/3) to linear partial equations of higher order and with 
constant coefficients and the characteristic features of 
their solutions have been extensively treated in the 
mathematical literature. We mention especially the 
work by Herglotz 38 and a recent investigation by 
Gfirding. 39 It seemed worth while to us to present a 
short account of some of Gardings results, since from 
these it follows that for all equations studied in III-A 
we have the same kind of relations between the initial 
value problem and the source problem as for (72), while 
moreover all these equations have propagation char¬ 
acter. In subsection 3 we will see that for equations of 
the exponential type (71) the situation is radically 
different. 


38 G. Herglotz, Ber. Sachs. Akad. Wiss, Math.-Phys. Klasse 
78, 41, 287 (1926); 80, 69 (1928). 

39 L. G&rding, “Linear partial hyperbolic differential equations 
with constant coefficients,” Acta Math, (to be published). See 
also Comptes Rendus 228, 731 (1949). It is a pleasure to thank 
Dr. Garding for making available to us a copy of his manuscript 
and for many instructive discussions on these topics. 
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2 . On the Propagation Character of Equations 
of Finite Order 

Girding considers polynomials <?(£i, • • •, £») in n 
variables & and with complex coefficients and associates 
with them partial differential equations of the form 

q(d/dx i, • • •, d/dxn)p=—p. (80) 

The first question is whether it is possible to determine 
solely from the structure of the polynomial the main 
characteristics of the solutions of (80). Especially one 
would like to know for what polynomials the equation 
is of hyperbolic character. Here it should be pointed 
out that an equation is of hyperbolic type only with 
respect to certain directions in the space £i, • • •, 
That is to say, only for surfaces perpendicular to these 
directions are the problems (a) and (0) correctly set 
and only with respect to these directions do the solutions 
have propagation character . 40 Thus for the wave 
equation these directions are the time-like vectors and 
correspondingly the surfaces on which initial conditions 
are specified are space-like surfaces. 

More precisely, one has therefore to find a criterion 
for a polynomial q(£ i, •••,£„) to be hyperbolic with 
respect to a fixed direction (|i, •••,{»). Gfirding’s 
criterion is: Put q—p+r, where p is the principal part 
of q, i.e., the homogeneous part of q of highest degree. 
Then if for a direction f* we have 

P( q(rii+irii)^ 0 , (81) 

for all r which are larger than a real number r 0 and for 
any real vector 77 *, then the polynomial g(£ t ) is said to 
be hyperbolic with respect to | t -. Gfirdings main theorem 
is: a necessary and sufficient condition for the solutions 
of (80) to have propagation character with respect to 
It is that qfe) is hyperbolic with respect to | t -. 41 

40 The other main property mentioned above, viz., the existence 
of a relation between the solutions of (a) and (/3), is well known 
[see e.g., Courant-Hilbert, Methoden der mathematischen Physik 
(Dover Publications, New York), Vol. II, p. 165] and is of a 
much more elementary nature. Suppose (80) to be of mth order 
in d/dt where the t axis is the hyperbolic direction. Con¬ 
sider then the initial value problem q<p=0 with initial con¬ 
ditions that for t= — oo all <p and all its time derivatives, 
up to the (m— 2)nd one are zero, while d m ~ 1 (p/dt m ~ 1 =p(s, x) 
where s has to be considered as a parameter. Let the solution be 
<p(t,x;s). Then the solution of the source problem for (80) is 
W, x) = f-J<p{t-s, x; s)ds. 

41 This is stated very loosely. More explicitly we mean by 
propagation character of (80) in the direction |* the fact that for 
planes perpendicular to this direction one can suitably generalize 
the problems (a) and {&) and that the solutions have again a 
finite domain of dependence. Garding defines propagation char¬ 
acter still more rigorously as follows: Take a sequence ^(zr • •£„) 
of solutions of (80) which is such that for k— tends 

strongly to zero in a plane lizH-+|nX„=0. (By “strongly’’ 

is meant that as well as all its derivatives tend uniformly to 
zero on every finite domain in the plane). If from the fact that 
\f/k tends strongly to zero in this plane it follows that «fa tends 
strongly to zero everywhere, then (80) is said to have propagation 
character with respect to j i. It is intuitively clear that if the 
solution of (80) has a finite domain of dependence the above- 
mentioned limit property is true. The converse is not so obvious 
and is shown by Girding by explicitly constructing the domain 
of dependence, to which question we will return presently. 


It is easy to apply (81) to the wave equation (72). 
In fact (72) can be considered as the equation adjoined 
to the polynomial 

q=z i 2 +£ 2 2 +£ 3 2 -£4 2 . 

taking |»= (0, 0, 0, 1); i.e., the direction of the time 
axis, one verifies (81) with t 0 =0. The same is true for 

q= £i 2 + & 2 + £ 3 2 ?4 2 - k*-»Q- k*)P= - p 

for real k. If k is complex the time axis is again a 
hyperbolic direction, but for r 0 we now have to take 

n = {i((a 2 +/3 2 )*—a)}*, K 2 =a+ip. 

The significance of r 0 will become clear presently when 
we will deal with the source problem (see reference 43). 

It may be remarked finally that also the polynomial 
to which we can adjoin 

no-«w —p ( 82 ) 

t 

is a hyperbolic polynomial for real or complex, single 
or multiple occurring This is true because (for 
suitably chosen r 0 ) each factor of the product is 
hyperbolic with respect to the time-axis. 

Consider now the source problem for (80). It is clear 
that a solution in terms of Fourier integrals can be 
given in the form 

1 /• /» dkdo) 

W, x ) = —— I I :—77 expp(kx-o^)] 

(2tt) 4 *7 J zk) 

xf (dt'dx r p(t\ x') exp[— i(kx'— a/)] (83) 


dt'dx'p(t\ x')K(t-t f , x-x 7 ), (84) 

1 r r dkdu 

X ) = 7TT: I I —;—— expp(kx-c^)], (850 

(27r) 4 J J q(—to3,tk) 

where in view of the physical applications, we consider 
from now on polynomials q in the four variables 42 x, t } 
which are hyperbolic with respect to the time direction, 
i.e. | t = (0, 0, 0, 1). It is also well known how to take 
care of the retarded action of the source by taking the 
path of the co-integration above the real axis in the 
co-plane such that all the zeros of q{— fco, ik), the latter 
considered as a polynomial in co, lie below this path. 
This recipe can also be expressed otherwise; viz., by 

42 We have followed the customary physical convention of 
taking opposite signs in the Fourier development in x and t. 
This is contrary to what is done in reference 39, and as a result 
certain differences in sign between our account and Garding’s 
paper appear. 

We have also with the usual lack of qualms interchanged 
orders of integration in going from (83) to (84). K(t, x) is the 
solution corresponding to a source 8(x) = 8(t) ■ 8(x). The careful 
justification of the result (84) is dealt with extensively in reference 
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putting 


K(t, x) = —fdk f du 

(2w) 4 J 


exppkx+ (r— io))t] 
q(T—iu, 7k) 


(85.) 


where r is an arbitrary quantity larger than some fixed 
real number r 0 and where we now integrate along the 
real w-axis. 43 Since K{t , x) is independent of the value 
of r, one sees from this form immediately that K(t, x) = 0 
if /<0 by letting r—©. 

The essential step in the derivation by Gfirding of 
the domain of dependence lies in a generalization of the 
form (852) to 


m ., x)= 


-i-r 

(27t) 4 ‘/_ c 


dk 


f 


X I do) 


exp[(r?+ ik)x+ (rjo- to)/] 

?(t£ o—to, r^+ik) 


, (85.) 


where again r>r 0 and where moreover a four-vector 
(Jo, ?) has been introduced which satisfies the require¬ 
ment that the zeros of q still lie below the real co-axis 
and all paths of integration are now taken along the 
real axis. From Cauchy’s theorem it is then clear that 
(85 3 ) is equivalent to (85 2 ) and that therefore the 
integral in (85 3 ) is independent of the vector & and of 
the magnitude of r. By letting again r-o one now 
sees that in all points t, x for which 

£o/+?x^O (86) 

we shall have K(t, x) = 0. To find therefore the domain 
of dependence of the solution of (80) one only has to 
determine the set of directions & which satisfy the 
condition stated above. For the equation (| | — n 2 )\p— p 
one readily finds, whether k is real or complex, that 

the Jrdomain is 

?o>0, to 2 -e>o, 

i.e., the inner part of the future light cone. Then it 
follows from (86) that the domain of dependence is 


part p of q and put 

^(rf+f)-^(f)n(r+«i({,{)). (87) 

X 

For real J, the roots Wi(J, J) can be shown to be real if q is hyper¬ 
bolic. From the homogeneity of p one easily sees that 

a£) = aui( f, J); u t (al J) = (1 /o)ui(£, £) (88) 

and 

«.4£) = 1. (89) 

The cone F(J) is defined by the assembly of all directions J for 
which the smallest of the roots m is still >0. From (88) it follows 
that we really have to do with a cone determined by the direction 
of J only, while from (89) and the continuity of J) one infers 
that such a cone must exist and that the direction J belongs to it. 
The importance of the cone T(J) lies in its following two main 
properties: (1) If q is hyperbolic with respect to J then it is also 
hyperbolic with respect to any direction J in F(J); and any 
direction in T(J) generates the same cone. (2) The cone F(J) 
contains all vectors for which (85 3 ) is identical with (85 2 ) and 
therefore according to (86), F(J) determines the domain of 
dependence. 

3. Remarks on the Propagation Problem for Equations 
of Infinite Order 

Little is known about equations of the type (71). 
However, it can readily be seen that there are essential 
differences from the hyperbolic differential equations 
discussed above. 

In the first place it is clear that in the homogeneous 
equation the exponential factor does not play any role 
so that the initial value problem is the same as for/=0. 
Thus with regard to the Cauchy problem the equation 
has propagation character. Another way of expressing 
this is to say that the D-i unction in (74) and (79) is 
unaffected by the introduction of an exponential. 

Just as for the exponential case in the mechanical 
problems (see II-B-3) the situation is quite different 
for the inhomogeneous problem. The solution of 

e /fa>Q/,= _p (90) 

can still be written down in a form analogous to (83): 

1 r r dkdoi 

'/'(/, x)=-I | --exp[—/(o) 2 —k 2 )+*(kx— 0 >/)] 

(2ir)0 J k 2 -to 2 


1^0, / 2 -x 2 >0. 

Clearly the same is true for a product of factors O - *%) 
so that we may state: The equation F([~M=P, where 
F is an arbitrary polynomial with complex coefficients, 
is a hyperbolic equation with respect to the time 
direction. Its domain of dependence is the past light 
cone and its interior. 

Garding has succeeded in expressing these results in a much more 
general way. He associates with the polynomial q which is hyper¬ 
bolic with regard to the direction j i a cone r(J t ) generated by the 
direction J t -. T(J) is defined as follows. Consider the principal 

43 Clearly, the meaning of this transformation is that since the 
number of zeros of q is finite we can push them all below the real 
w-axis by adding an imaginary part to w; the quantity ro which 
we have encountered in Eq. (81) and in the examples discussed 
above is simply the imaginary part of that root of q(~ia> f /k)== 0 
which lies highest above the real co-axis. 


X J Jdt'dx'p(t ', x') exp[—7(kx'~coO] 


(91) 


where one would be inclined, as in subsection 2, to 
perform the oo-integration along a path parallel to 
and above the real axis. However, the argument that 
this choice will give a retarded action breaks down since 
one clearly may never close the path by a big half 
circle, so that one cannot conclude that \p=0 for £<0. 

Nevertheless it is of interest to discuss the formal 
solution (91) for this choice of path. For the case that 
/O) is a polynomial one still must make an important 
distinction, viz., that between odd and even functions 
/. This will be clear enough from the two cases we will 
now discuss separately: 

(a)/0)=-X 2 Q (b)/G)=XO, 
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where X is a quantity with the dimensions of a length. 
The signs have been chosen such that for the static 
point source problem, (90) always gives a \p finite at 
the origin. 

Case a .—If p is the 5-function 5(/)5(x), the w-integral 
is divergent so that the analog to D ret (x) does not exist. 
This does not mean that (91) always leads to a divergent 
\p, but it implies a severe restriction on the admissible 
source functions p. This can also be illustrated as fol¬ 
lows. For the case considered (90) becomes 

exp(—XO -QKi, x, X 2 ) = — p(t, x), (92) 

where the dependence of ^ on the parameter X has been 
written down explicitly. Hence differentiating with 
respect to X 2 and omitting the plane wave solutions for 
i/' which do not interest us at the moment, we get 


Now divide the ( k , gj)- domain in two regions and make 
the following substitutions: 

k 2 >oo 2 : k — R coship, w=jRsinh<p; — <x><R<+ cc 
k 2 <cc 2 : k—Rsmhtp, u = R coship; — oo<<p<+oo. 

One then obtains after some calculation 


1 1 \R\ 

7=- - dR exp(-X 4 i? 4 ) d§ 

47T 3 R F J 0 

X {sinh# sin(pbR sinhi?) — cosh$ sm(plR cosh#)} 


1 d r™ dR r w 

= -—exp(—X 4 jR 4 ) I d# 

7 T 3 dpJ o R J0 


X {cos(pbR sinh#) — cos(pbR cosh??)}, 


#/a(X 2 )-O=0. (93) or finally 44 


Equation (93) has to be solved with the initial value 
\p(t, x, 0) satisfying the ordinary wave equation 

OKt, x, 0) = — p(t, x). 

Equation (93) is of the heat conduction type in the 
variables X 2 and x but with respect to X 2 and t it is the 
analog of the equation of heat conduction toward 
negative time. Now it is well known that for the 
ordinary heat conduction equation 

du/dt=d 2 u/dx 2 

the solution is an analytic function of x for any />0, 
whatever function the given initial value u(x, 0) is. 
Clearly therefore, in order to get a solution for an 
equation of the type 

du/dt— —d 2 u/d0i c®, 

one can only allow initial values which are analytic 
functions of x , and even this may not be enough. 

Case b .—Here (91) exists even for a 5-function source. 
To evaluate it further one conveniently cuts up the 
path of the co-integration into: (1) two small half-circles 
over the poles co=rb|&|, (2) a principal-value integral 
along the real axis.- One easily sees that, as in the case 
X = 0, the contribution of the half-circles is 

-(l/8irr){5(/+r)-5(/-r)}. (94) 

The principal value integral I becomes, after some 
angular integrations, 

1 1 r* r™ exp[—X 4 (& 2 —co 2 ) 2 ] 

7=-I kdk-P I dco - 

16iw 3 rJ _oo k 2 —co 2 

X [expp ( kr —co/) ]—exp[ — i {kr+ co/) ]]. 

To evaluate this further, one has to distinguish several 
cases. Consider first />0, t 2 >r 2 and put 

/=p*cosha, r=p^sinha, p=/ 2 —r 2 . (95) 


1 d dy f ir 1 

7=-- exp(—X 4 y 4 p~ 2 ) K 0 (y)+-Y 0 (y) . 

7r 3 dpJo y L 2 J 

One proceeds likewise in the other parts of the r, /-plane. 
Adding the contribution (94), the final result for is: 

1 

<A=-[$(<+>■) —S(<—c)]+w(p) (96) 

87rr 

1 d l r™ dy 

<P\(p)= -1 sign(p) • I — exp(-Xyp~ 2 ) 

TT 3 dpi y 

x|*o(y)+^Fo(y)]j. (97) 

For X=0, ip can be written as 44 

,°(P) = ^(^ sign(p)^ • JT ^iSr 0 (y)+^F 0 (y)J 

2 z* 00 dy r* 3 du y 2 

= —5(p) I — I —sinw sin— 

tt 3 Jo y Jo u 4u 

1 r™ du dv v 1 

= —5(p) I —sinw I —sin-=—5(p). 

tt 3 J 0 u Jo v u 4 t 

Therefore, remembering the definition (95) of p, we 
have verified that for X=0, ^ becomes the function 
D iet of (76). 

For Xt^O, <px(p) is zero for p=0 and has two peaks 
symmetrically around p=0 at a distance ^X 2 , with 
height ~X~ 2 and width ^X 2 . For decreasing X the peaks 
will move nearer together and become higher, but will 
remain separated by the zero value of ip at the origin, 


44 K 0 and Fo are defined as in Watson, Theory of Bessel Functions 
(1944), second edition, Eq. (4), p. 78, and Eq. (2), p. 64, respec¬ 
tively; see also Eq. (14), p. 183, Eq. (13), p. 180, and Eq. (3), 
p. 184. 
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however small X may be. Therefore the limit at which 
the two peaks coalesce to the 5-function is reached in a 
highly non-uniform manner. 

Thus according to (96) the \p corresponding to a pulse 
(5-source) at the time /=0 is for />0 not strictly zero 
outside the sphere r=t, although for p>\ the excitation 
function \p drops rapidly to zero. Moreover it follows 
from (96) that for times prior to the occurrence of the 
pulse there is already an excitation. For all these 
reasons the notion of propagation therefore now loses 
its meaning. 

If, however, we average yp over a region of the order 
X around the light cone this average value is vanishingly 
small for />0, while for /<0 one approaches rapidly 
the same value D ret for yp as in the conventional theory. 
In this sense we can therefore say that the equation 

exp(X 4 Q 2 ) • O = - p (98) 

exhibits a propagation character in the mean . 

From the examples discussed it is clear that for the 
exponential type of Eqs. (90)—or, for that matter, for 
Eq. (71)—the propagation character in the sense dis¬ 
cussed in subsection 2 does not seem to exist. On our 
request Dr. Girding has kindly investigated this prob¬ 
lem along the lines mentioned in reference 41 and has 
given a mathematically rigorous proof of the non¬ 
existence of the propagation character for the Eqs. (71). 

For the purposes of the next section we note that it 
is still possible to define a function D\ for Eq. (98) in 
much the way that was done in (77) for X=0. One can, 
in fact, determine the solution of the equation 

exp(XO) -QH - 5(x) (99) 

which corresponds to the advanced solution for X = 0 
by integrating below instead of above the poles on the 
real co-axis. Defining then D\(x) as half the sum of this 
solution and of (96), one obtains 

D x (x) = <p\(p). (100) 


Finally one can integrate (99) below the pole at 
03 — | k | and above the pole co=+ | k |. This mode of 

integration has been used extensively by Feynman in 
his version of quantum electrodynamics. For this 
integration the result is 


or 




i 

- 

4tt 2 


1 


r 2 —t? 


P~D x (x)+iiDV(x), 


( 101 ) 


where D a) (x) is the same function as occurs in the 
customary formulation of quantum electrodynamics. 45 
Mathematically speaking, the reason that D (l) comes 
out unchanged is that it is a contribution from the 
poles at co==dh |k| and is not affected by the behavior 
of \p at the point at infinity in the co-plane. 


46 See reference 37, Eq. (A.40). 


4. On Exponential Modifications of Quantum 
Electrodynamics 

In this section we shall discuss some of the novel 
features which arise if we modify the electromagnetic 
field equations in the manner indicated by Eq. (90), 
so that we start from a Lagrangian density 

1 dA M dA M 

£== - e /(D) - j»A,i—'t'(y»dn+K)'h (102) 

2 dx v dx v 

where all quantities are defined as in Eq. (52). We will 
again impose the supplementary condition (53) on A M . 
As £ exhibits all conventional invariance properties it 
is clear that an energy-momentum tensor density must 
exist which satisfies the usual conservation laws. (The 
most suitable way to find this tensor is by means of an 
elegant method recently developed by Green. 46 ) 

We consider first the case in which there is no 
interaction. It follows directly from the results in 
subsection 3 on the initial value problem that classically 
the situation is unchanged; i.e., it is the same as for 
/= 0. This is true also quantum mechanically. 

First of all we can make a spatial Fourier decompo¬ 
sition of the electromagnetic field quantities, as a result 
of which the electromagnetic field Lagrangian becomes 
a superposition of harmonic oscillator Lagrangians 
each modified by an extra exponential factor 
exp[f(—D 2 — co 2 )], where 03 is the frequency of the wave 
considered and D again stands for d/dt. From the 
argument given in II-B-3 it then follows that the result 
of quantizing each Fourier component is in all respects 
identical with that of an oscillator without exponential 
factor. Hence it follows also that the commutation 
relations between two field quantities taken at the 
same time but at different points in space remain 
unchanged. But, as according to subsection 3 the 
propagation character for the homogeneous problem is 
unaffected by the presence of exponentials, the commu¬ 
tation relations are also unchanged if we take the field 
variables at different times, so that we may state that 

[^4 m (x), A v {x')~\ = ib iiV D(x— x'), 

where D is the D-function of ordinary quantum 
electrodynamics. 

As further consequences of the absence of any 
influence of the exponential on problems in which no 
interactions are involved, we mention that Planck's 
radiation law is unaffected and also that the vacuum 
expectation value of any function of the A M remains 
unchanged. In particular, 

({ A M (x), A v (x ')} ) vac = 5 Aty Z> (1) (x—x') 

in which the conventional Z) (1) -function occurs. This is 
in harmony with the result quoted in Eq. (101) where 
in other context we met an unmodified Z> (1) . 


46 H. S. Green, Proc. Roy. Soc. 197, 73 (1949), Section 4. 
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Turning now to problems of interaction, a remark on 
the classical theory is in order. Calling <p the modified 
Coulomb potential, it is seen from (102) that for a 
point electron at rest 

e r™ dk 

<p(r) = -I — sinkr exp[/(— k 2 )~]. 

k 


Kallen and by Yang and Feldman. 47 Their idea is to stay 
in the Heisenberg representation throughout the calcu¬ 
lations and in this representation to solve by iteration 
the equations describing the behavior of sources and 
fields. 

We consider first the electron self-energy integral 
which in the notation of Schwinger is given by 48 


For both cases a and b discussed in subsection 3, the 
^-integral is convergent and tp is finite for r— 0. There¬ 
fore the classical self-energy of a point charge is finite, 
and this is true in any Lorentz frame of reference. For 
/= —X 2 Q this has been discussed by Born. 19 Thus for 
such static and quasi-static effects there is no essential 
difference between the cases a and b On the other hand, 
the truly dynamical features have been shown there 
to be quite different, and for the reasons given in 
that section we will from now on confine our atten¬ 
tion to the case of even functions / and will choose 
again especially /= X 4 0*. 

For this case we have given in subsection 3 the 
formal solution of the pulse problem, i.e., the excitation 
which is generated by a 6-function source acting at 
/=0, x=0 has been computed. It was found that such 
pulses act in such a peculiar way, viz., “before 
they start”, and giving an excitation which partly 
moves with velocities greater than that of light, that it 
is out of the question to consider without further ado 
the application to actual physical situations. In partic¬ 
ular the integrability conditions of the Tomonaga- 
Schwinger theory which express the unambiguous 
causal development with time of the state vector of a 
system can no longer be expected to be satisfied. 

On the other hand, we have been stimulated by the 
reasonable behavior in an average sense explained in 
the previous section to make some conjectures as to the 
possible use of such exponential equations in a situation 
where one no longer insists on the unlimited localiza- 
bility of space-time events. The tentative and incom¬ 
plete nature of the remarks made hereafter need, of 
course, hardly be stressed. 

First of all, it is clear that the effect of the modifica¬ 
tions may be quite different for virtual as compared to 
real processes. In the latter case the lack of propagation 
character is such a drastic change of our fundamental 
physical notions that one can only hope to obtain an 
interpretable scheme by renouncing a strictly causal 
description. This will be the subject of further investi¬ 
gation. 

On the other hand, one may expect that the effect on 
virtual processes will simply be the replacement of the 
ordinary D- and A-f unctions by the corresponding 
exponential modifications of these propagation func¬ 
tions, while the Z) (1) - and A (1) -functions remain un¬ 
changed. In fact, this procedure can be justified by 
means of the 5-matrix methods recently developed by 




+5(£)D <1) ($)]y„ exp 


In keeping with the foregoing, one will retain the 
Z) (1) -function but must now replace D by D\ given by 
(97) and (100). This would make the first term in the 
square brackets finite, but would not alter the second 
which still diverges. Thus, exponential modifications of 
the electromagnetic field do not yield a finite self-energy 
in contrast to the polynomial modifications studied in 
previous sections. This is essentially because the expo¬ 
nentials do not affect the vacuum fluctuations, whereas 
in the polynomial case new fluctuations occur due to 
the presence of new types of quanta. 

It would seem plausible, however, that once one gives 
up the localizability of the electromagnetic field quan¬ 
tities, the same should be done with the electron field. 
Thus, one will keep again the expression for A (1) but 
may let A satisfy 

exp[X 4 0- k 2 ) 2 ] • O-k 2 )A(*) =-S(x), (103) 

where the exponent has been so chosen that according 
to a similar argument as given above the free electron 
states are unaffected. This implies a modification of 
the last term of (102) by the exponential factor 
expX 4 (|~~]— * 2 ) 2 . Thus, the matter equations now are 

exp[X 4 0— k 2 ) 2 ] • (y A+ (104) 

Doing this, one obtains an absolutely convergent 
integral for 8 m which can be readily seen to be 

€ 2 /c z* 00 r 1 [ 2—(l+y 2 )* 

8 m =— I y 2 dy\ - \ - 

87tV 0 La+y^M-l+a+y 2 )* 

2-f- (1-f-y 2 )* 

Xexp[—4X 4 k 4 (1 — (1+y 2 )*) 2 ]- 

l+U+y 2 )* 

I 2 

Xexp[—4X 4 k 4 (1+ (1+y 2 )*) 2 ] H— exp(— 4X 4 * 4 y 2 ) 

) y ' 

Once the A-function has been modified according to 
(103) the question arises as to how this will affect the 


47 G. Kalian, “Mass and charge renormalizations in quantum 
electrodynamics without use of the interaction representation 
Arkiv. f. Mat. Astr. o. Fys. (to be published); C. N. Yang and 
D. Feldman (to be published). We are indebted to Professor 
Pauli for making Dr. Kallen’s manuscript available to us as well 
as for instructive discussions on this point. 

48 See reference 37, Eq. (3.78). 
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vacuum polarization. The integral representing the 
current h{ie$y^) induced by an external field A / is, 
according to Schwinger, 49 

d(ie\py^l/)= — 4:e 2 JdiX f K tltf (x—x , )A v € (x , )j (105) 

dA dA (1) dA dA (l) 

K^(x) = -1- 

dXn dx v dx v dx» 

/ dA dA a) \ 

-Sj - +k 2 AA^ ). 

\ dx\ dx\ / 

If instead of the ordinary # A-function one again employs 
the expression following from (103) it is not surprising 
that the charge renormalization will be finite. The more 
interesting question is what will now happen with gauge 
invariance and charge conservation. The decisive quan¬ 
tity here is the divergence of for which one has 

dK^/dx^ (dA^/dx,) • (□- k 2 ) A 

= — (dA a) /dx v ) *exp[—A 4 0— k 2 ) 2 ]$(x) . (106) 

This is not equal to zero as one would like it to be. 
Nor could this be expected since (104) is not gauge 
invariant. This of course makes it questionable to call 
a current and the quantity (105) an induced 
current. It can be asked, however, whether (106) 
could be zero “in the mean”; i.e., in some way similar 
to the manner in which we obtained an average propa¬ 
gation of an electromagnetic excitation with light 
velocity. For this purpose we study 

x(*) = exp[-X 4 0- <c*) 2 ]s(*). 

By the same methods used to derive Eq. (97) one finds 


1 d 1 r“ - 

r // Vi 

x(*) = -I ydy ex 

p -X 4 ( -+K 2 ) . 

7T 3 dp pJ o 

L \ p / J 


We note the following main properties of x: First it is 
when integrated over a space-time domain ^X 4 
around the origin and in this sense behaves approxi¬ 
mately like a 5-function. Second it vanishes very 
strongly for /=0, x=0. In fact, the product of xM and 
the derivative of A (1) occurring in (106) is zero at this 
point, notwithstanding the singular behavior of A (1) for 
p=0. Thus one can give an unambiguous meaning to 
the right-hand side of (106), in contradistinction to the 
case that X=0. As x(%) is an even function of its 
argument while dA a) /dx v is odd, it follows that the 
average of (106) over a small space-time domain around 
the origin vanishes. In this sense we can speak of a 
conservation of the induced charge in the mean. 

49 See reference 37, Eq. (2.20). 


The preceding remarks would seem to indicate that 
the assumption of non-localizability should be made for 
all kinds of fields once it is made for one kind, and that 
“macroscopic” conservation laws may perhaps have to 
be considered as valid only in some average sense. 

IV. OUTLOOK 

From the investigations described in this paper we 
have seen that it is very difficult, if it is at all feasible, 
to reconcile the three requirements of convergence, 
positive-definiteness and strict causality in dealing with 
field equations which are partial differential equations 
of high order or are integral equations. It may be asked 
whether the restriction to integral functions made in 
this paper has been unduly narrow and whether more 
general functions might be better suited for the purposes 
outlined in the introduction. We have not made a 
systematic exploration of such possibilities. A further 
study may be of interest, but it seems to us that the 
main physical consequences of the use of finite distance 
operators have been laid bare by the present investi¬ 
gation. However one point which we still wish to 
mention is that by using meromorphic functions or 
rational fractions for F in (2) it is possible to obtain 
multi-mass equations without having negative-energy 
troubles. From Eq. (47) the reader will in fact verify 
that if between each two zeros of F there lies a simple 
pole, one gets an assembly of fields with quanta of 
various masses and all with positive-definite energy. 
On the other hand, the introduction of fractions leads 
again to loss of convergence. 50 

The multi-mass equations studied in this paper are, 
generally speaking, connected with the idea of compen¬ 
sation (see Introduction). As we have seen, the various 
fields whose contribution to the self-energy etc. tend 
to balance each other are all of the same relativistic 
transformation type. Other compensation schemes (see 
reference 6), are not covered by the present work. On 
the whole it seems however that all attempts made so 
far to achieve convergence by compensation are inade¬ 
quate. This is in our opinion not so much a criticism of 
the underlying idea as well an indication that one should 
look for still more intimate connections between fields 
of various types. 

All in all, it seems to us that the most hopeful 
approach which the present investigation possibly indi¬ 
cates is one, where, as for example in the case of the 
exponentials studied here, the possibility of an ordering 
of space-time events is no longer a strict requirement. 
But it should again be emphasized that here we have 

50 To quote an example: consider 

□ * tan\D*-^ = F(D)^==-X; M . 

For X—►O (or for low frequencies in the Fourier development) 
F(m)^XD. F has zeros k„ 2 , where —wirX” 1 , n = 0, 1, 2, ••• 
are quantum rest masses. Also F has singularities at □=a, l 2 , 
«„=*(»—I)irX -1 . Thus, all quanta have positive definite energies. 
Convergence is lost as is already clear from the singular nature 
of the modified Coulomb potential which is (47ir) -1 coth7ir/2X and 
is ^r~ 2 for r—*0. 
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come to the limits of the foundations of our present 
picture of the physical world. An attempt at loosening 
up these foundations by relinquishing the “causality 
in the small” might throw more light on the situation. 

The present work was begun during the stay of one 
of the authors (G.E.U.) at the Institute for Advanced 
Study. He wishes to express his appreciation for this 
stimulating experience and for the hospitality shown to 
him. He is also indebted to the Institute for having 
supported his research in part by a grant of AEC funds. 
Both authors also wish to thank the staff of the Brook- 
haven National Laboratory for working facilities pro¬ 
vided them last summer. 

APPENDIX 

Some Further Comments on Multiple Frequencies 

We wish to indicate here briefly the treatment of a mixture of 
single and multiple frequencies, and that of frequencies with 
multiplicity higher than two. For this purpose we consider first 
the situation in which three distinct real frequencies are present: 


2 


z.= -? n 

£** 0 

(Al) 

^rnQiilP+uftQr, 

0 

(A2) 

i 


770 = {(cdl 2 — W0 2 )(cO2 2 — COO 2 ) I -1 , 

7/1= {(o>2 2 — «l 2 )(wi 2 — COO 2 ) }“ l , 

772 = { (cx> 2 2 — Wl 2 )(cu 2 2 — COo 2 ) l“ l . 



From this we get the case of one single plus one double root by 
putting co 2 =«i+e } Qi = (> 2 -f*e (?3 and going to the limit e-*0 in the 
same way as described in II-A-4. Obviously the Co-dependent 


part of (A2) remains unaffected apart from i 7 o-K«i 2 —"o 2 )~ 2 . The 
generalization to the case of more than one simple frequency plus 
multiple frequencies is obvious. We can also get from (A2) the 
case of a triple frequency by putting 

Wl—WoT^l, £02 = COo-f-€2, 

Co = ^0+(ei + €2)Rl + K c l4-€2) 2 ^2, 

Qt — Pod - eiRi+lei 2 /^. 

The transformed L then becomes (after multiplication with a 
trivial constant factor — 4w 0 2 ) 

L — RqR% — COo 2 Ro2?2 — Rl 2 — Wo 2 Rl 2 T4<jdoi?0^l} 

which can be shown to give the same equations of motion in 
terms of q as (Al) does. The Hamiltonian is 

H = PoPz-^iixpRaRz-h 4Pi 2 -f-coo 2 Pi. 2 4-Po 2 —4co 0 RoRi. 

Performing the contact transformation 

Po—cooHPo- 7Z rrPiy, Ro~~* ;Po 

V coo(Aoo)* / coo* 

P^fp.+^-P.), Rl ^ Rl+ -±-^R^ 

P 2 — -R, 

\ ( 2 coo)* / coo’ 

// becomes 

H**U»(P 0 R,-P t Ro)+-^-W + Rt)+\(Pl t +a> l ?R l *). 

4 coo l 

Comparison with (37) shows that a triple frequency is equivalent 
to a double one plus a single oscillator with the same frequency 
which of course introduces a new degeneracy. 

The generalization of relations like (A3) to the case of w-fold 
multiplicity is straightforward. We have not investigated this 
further but conjecture that a 2w-fold root is equivalent to a 
degenerate system of n double roots and that a (2»+l)-fold root 
is equivalent to a 2«-fold root degenerate with a single oscillator. 



